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ABSTRACT 



Context. In the standard formulation of the disc instability model for soft X-ray transient and dwarf nova outbursts, the mass transfer rate 
from the secondary is assumed to be constant. This may seem natural since the L\ point remains shielded from the accretion luminosity by 
the accretion disc. However an indirect heating could take place and could lead to an enhancement of the mass transfer rate. It is still debated 
whether such an enhancement of the mass transfer rate during outbursts is a missing ingredient of the current model. 

Aims. We discuss two mechanisms that could result in an enhancement of the mass transfer rate during outbursts: the hot outer disc rim itself 
could significantly heat the L\ point; scattered radiation by optically thin outflowing matter could also heat L\ significantly. We determine 
quantitatively the increase of the mass transfer rate resulting from an extra heating. 

Methods. We estimate the heating of the L\ point by the disc rim by analytical and numerical arguments. The fraction of the luminosity 
scattered by matter above the disc, here modeled as a spherically symmetric outflow, is determined analytically. We finally solve in a numerical 
way the vertical structure equations for the secondary star and calculate the mass transfer enhancement. 

Results. During outbursts, the temperature at the outer edge of the disc reaches 10 4 K. The disc edge heats up the upper layer of the secondary 
with a flux of the order of the intrinsic stellar flux. This probably has no large effect on the mass transfer rate. In soft X-ray transients, the 
environing medium of the disc (corona+wind) could back-scatter a certain fraction of the accretion luminosity toward L\. In dwarf novae, the 
same effect could be due to the wind present during outburst. Since soft X-ray transients reach high luminosities, even a low efficiency of this 
effect could yield a significant heating of L\, whereas we show that in dwarf novae this effect is negligible. Initially the incoming radiation 
does not penetrate below the photosphere of the secondary. Depending on the heating efficiency, which has to be determined, the mass transfer 
rate could be significantly increased. 
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1. Introduction 

Dwarf novae (DNe) and Soft X-ray transients (SXTs) are semi- 
detached binary systems that undergo regular outbursts. SXTs 
in outbursts reach luminosities of order of 10 38 erg.s -1 in the 
soft X-ray band, whereas they are faint X-ray sources during 
quiescence (with luminosities < 10 33 erg.s -1 ). Their outbursts 
last for a few mon ths and have a re currence time in the range 
1-50 years (see lLewin et al.lll997l) . DNe are 1-5 magnitude 
brighter during outbursts than in quiescence. Their outbursts 
last f or a few days and recur every few weeks (see IWarneii 
2003b . In both types of systems, a compact object (i.e. a white 
dwarf in DNe and a neutrons star/black hole in SXTs) accretes 
from a low mass secondary star via an accretion disc. It is 
now believed that DNe and SXTs outbursts are both due to a 
thermal/viscous instability of the accretion disc triggered when 



hydro gen becomes partially ionized (see JOsakil 1 1 996c lLasota 
200 lh . In this picture, the accretion disc performs a limit cy- 
cle between a cold, quiescent state of low accretion rate and a 
hot, viscous state of high accretion rate corresponding to the 
outburst. The difference between DNe and SXTs timescales 
is thought to be due both to the irradiation of the disc by the 



centra l source in SXTs (see iKing & Ritten ll998: Dub us et al 
200 II) . and to a higher primary mass. 
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The thermal/viscous instability alone is not sufficient to 
explain the rich and di verse feature s of DNs and SXTs 
lightcurves. For example. lOsakil dl989l) su ggested that a tidal 
instability driving the disc eccentric (see IWhitehursj ll988) 
would account for the superoutburst phenomenon of SU UMa 
stars (a subclass of DNe) when coupled to the standard ther- 
mal/viscous instability. Tidal effects have also been proposed 
to account for t he secondary max imum often present in SXTs 
lightcurves (see lTruss et al.ll2002l) . 
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Table 1. Physical parameters for the four models condidered in this paper. i\r 1S the orbital period in hours, a the orbital 
separation, Ro the secondary radius, M\ and M^ the primary and secondary masses, T* the effective temperature of the secondary, 
M tI the average mass transfer rate from th esecondary , and M accmm i s the maximum accretion rate onto the compact object. T+ 



for model 1 - 3 is taken from Smak (2004) and from Ham eury et alj dl986l) for model 4 



Model 


System 


Port, (h) 


a (Re) 


R 2 /a 


q = Mi/ Mi 


M 2 (M e ) 


r*(K) 


M tl 


(10 16 g.s- 1 ) 


M acc , max (10 16 g.s- 1 ) 


h/r 


1 


OYCar 


1.51 


0.6 


0.21 


0.1 


0.085 


2500 




0.5 




16 


0.06 


2 


UGem 


4.24 


1.48 


0.29 


0.36 


0.7 


3500 




5,10 




100 


0.09 


3 


Z Cam 


6.96 


2.17 


0.35 


0.6 


0.7 


4200 




30,60 




250 


0.1 


4 


A 0620 


7.8 


3.9 


0.21 


0.1 


0.7 


3000 




0.5 




10 J 


0.1 



In its most common formulation, the disc instability model 
(DIM) assumes that the mass transfer rate from the secondary 
is fixed, and given by the secular losses of angular momentum 
from the binary system. This assumption is most likely to be 
wrong and significant deviation of the mass transfer rate from 
its secular mean on time scales as short as the outburst dura- 
tion could be an essential ingredient missing in the standard 
model. Numerical investigations showed that variations of the 
mass transfer rate could explain the outburst bimodality of DNe 
dSmakll 1999b a nd that it co uld lead to long outbursts similar to 
superoutbu rsts (Hameury 2000). Auguste iin et al.l y_993) and 
Chen et al.l d 19931) sought to explain the rebrightenings of SXTs 
lightcurves by an increase of the mass transfer rate. However, in 
these works either an empirical relationship between the mass 
transfer rate and the mass accretion rate was assumed or mass 
transfer burst episodes were imposed without physical grounds. 
It has been often argued that a such mass transfer enhancement 
could result from irradiation of the secondary during an out- 
burst. IVialleX&TIarneur^l (J2007|) investigated this problem and 
found that, because the L\ point is shielded from direct irradia- 
tion disc by the secondary star, the increase of the mass transfer 
rate is modest. 



secondary. The resulting heating is probably too small to have 
any major effect. For soft X-ray transients, the backscattered 
radiation is in the soft X-ray band and is therefore strongly ab- 
sorbed by the initially neutral hydrogen in the very upper part 
of the atmosphere. However, this large X-ray flux significantly 
affects the atmosphere. A full analysis of the radiative transfer 
should be undertaken to determine the exact structure of the 
atmosphere. 

The paper is organized as follow: in Section 2, we discuss 
the two mechanisms mentioned above. Section 3 is devoted to 
the computation of the vertical structure of secondary stars in 
DNe and SXTs, and to the relationship between the mass trans- 
fer enhancement and the heating of L\ . Concluding remarks are 
given in section 4. 

2. Heating flux at the L\ point during outburst 

In the following, we considere 4 models, with parameters given 
in Table 1, that are representative of various DNe subclasses 
and SXT outbursts. Models 1-3 correspond to dwarf novae with 
increasing orbital periods, and model 4 corresponds to A 0620- 
00, emblematic of short period soft X-ray transients. 



In this paper, we investigate two different mechanisms that 
could be responsible for a mass transfer enhancement, result- 
ing in both cases from an increase of the temperature at the L\ 
point. Since L\ lies in the shadow cast by the accretion disc, 
any heating must be indirect. The first effect that we consider 
is heating by radiation emitted by the hot edge of the disc dur- 
ing an outburst, which produces a flux at L\ comparable to the 
intrinsic stellar flux. The second effect that we consider is the 
heating of L\ by the accretion luminosity scattered by a wind 
and/or a corona extending above the accretion disc. The scat- 
tered fraction is small, but the scattered flux is still significant 
in SXTs due to their high luminosity in outburst. We then turn 
to the determination of the mass transfer enhancement. First, 
we compute the vertical structure of the secondary envelope in 
DNe and SXTs both in quiescence and outbursts to determine 
the dependence of the mass transfer rate versus a given incom- 
ing heating flux. Our results show that M can potentially be 
raised by a factor of up to ~ 100 in SXTs, and up to ~ 10 in 
DNe. We then discuss if the effects investigated here could lead 
to an efficient heating of the L\ point. This is done by com- 
puting the fraction of the incoming radiation that penetrates 
the photosphere. This depends on the nature of the radiation, 
hence on its origin. We show that the thermal flux emitted by 
the edge of the disc can penetrate below the photosphere of the 



2.1. Heating by the disc rim 

During an outburst, the disc grows as a consequence of an 
effective outward transport of angular momentum. The max- 
imum disc size is determined by tidal interactions with the 
secondary which effectiv ely trunca t e the disc at a mean ra - 
dius R tid ~ 0.S5R Ll (see lSmaril2003llchikawa& Osakilll994l) . 
During an outburst, the disc will reach this maximal radius, 
for SU UMa systems this would be true only during super- 
outbursts ("see IOsakilll989L Il996[). and tidal dissipation heats 
up the outer region of the disc. ISmakl (120021) argued that tidal 
dissipation occurs in a very narrow region at the disc rim and 
suggested that this energy is radiated away by the disc edge, 
whereas heat generated by viscous dissipation is radiated by 
the surface of the disc; this contrasts with our previous work, 
where tidal dissipation was considered to occur on a relatively 
large are a of the disc, significantly larger than the disc fhick - 
ness (see iBuat-Menard et alJuOOU lHameurv & Lasotal 120051) . 
We do not discuss this issue further here and we use Smak's 
assumption which maximizes the energy radiated by the disc 
edge. We show below that, even in this very favorable hypoth- 
esis, the effect is small. 

In this picture, removal of angular momentum by the tidal 
torque is negligible inside the truncation radius. The only effect 
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of the tidal torque is to prevent the disc from exceeding the 
radius R^. In this prescription for tidal effects, the only free 
parameter is R t ^. 

In the following, we determine whether h\ can be signifi- 
cantly heated by the hot edge of the disc. We consider this effect 
only for dwarf novae, since for soft X-ray transients the major 
effect in heating L\ is due to scattering by a ho t corona or a n 
outflowing wind (see the next subsection). As in lSmakl (12002b . 
we first determine the tidal torque Aid applied to the disc by the 
secondary, equal to the rate at which angular momentum is re- 
moved from the disc. The work done by this torque is radiated 
as a "tidal luminosity" L t a - 7 t id^, where Q is the Keplerian 
angular velocity at the disc edge. L t ;d is assumed to be radi- 
ated by the disc edge, whose effective temperature r e d ge is then 
given by: 



H\ 



Lad = <TTtj2xl&- 



edge 



r \r„ 



(i) 



where Rj is the disc outer radius and the ratio H/R at the outer 
radius appears explicitly. H is the full geometrical disc thick- 
nes s, of the orde r of a few times the vertical scale height. 



Smak (2002) considered the steady state case, which made 
the determination of Aid easy. As we are interested in Aid dur- 
ing an outburst, we consider instead the local balance of an- 
gular momentum at the disc edge: once the disc has reached 
its maximal radius, the angular momentum transported at the 
disc edge has to be fully removed by the tidal torque. Using the 
standard conservation laws for angular momentum and mass, 
this yields: 



Aid = (3ttvZ - M a )j d + j s M lr 



(2) 



where ja and j s are respectively the specific angular momentum 
at the outer disc edge and the specific angular momentum car- 
ried by the stream. The first term on the right hand side includes 
the viscous torque at the outer radius, i.e. the rate at which an- 
gular momentum is transported outwards by viscosity, and the 
effect of the incoming strea m (M tr is the mass tran sfer rate). We 
use the numerical code of lHameurvet al.1 dl998l) to calculate 
the time evolution of the accretion disc and thus determine Aid- 
The code has been modified to account for the different pre- 
scription for the outer disc radius. When R < R t g, we use the 
free outer radius conditions (Eqs. 7, 8 in lHameurv et al.l ll998, 
note that Eq. 8 is equivalent to .Aid = 0). When R reaches R t a, 
the code switches to a fixed outer radius condition R - R tl d- Eq. 
(f2| is then nonzero and yields the rate of angular momentum re- 
moval by the tidal torque without the need for tidal terms in the 
equations. The code switches back to the free outer radius con- 
ditions when 7 t id falls below an arbitrarily small, given value. 
We do not include here any heating of the outer disc region due 
to the impact of the incoming stream (see discussion in lSmakl 
2002h . 

Figure 1 shows the results for models 1 (left panel) and 2 
(right panel). In our model 1, the disc always fills its tidal ra- 
dius. Outbursts are of the inside-out type and heating fronts 
propagate over only ~ 80% of the disc. The heating front then 
dies out in the outer part of the disc and the corresponding out- 
ward transport of angular momentum is large enough to coun- 
teract the effect of the incoming stream. This contrasts with the 



Table 2. Results for the heating of the secondary by the edge 
of the disc. r e d ge ,max is the maximum effective temperature of 
the disc edge. F is the heating flux impinging the L\ point (F* 
is the intrinsic stellar flux). 



Model 



' edgejnax 



(K) F/F* 



9000 2 

13000 5 

10000 2.5 



norm al tidal torque prescrip tion model (as described for exam- 
ple in lHameurv et a l. 1998), where the outer radius of the disc 
is found to vary by about 10%; the reason is that the disc is too 
large - despite the short orbital period - given the strength of 
the outbursts that cannot reach the outermost parts of the disc. 
In Fig. 1, one can also see that the edge of the disc is heated 
for a longer time than the outburst duration. The effective tem- 
perature of disc edge is maximum during outburst and reaches 
10 4 K. For model 2, the disc outer radius varies, and the disc 
reaches its tidal radius during outburst. The effective tempera- 
ture of the edge reaches a maximum of 1.3 10 4 K during the 
outburst. Note that in order to obtain a sequence of regular out- 
bursts for this model, we used a mass transfer rate of 10 17 g.s 
~' ; for a value half the size, we obtain an irregular sequence of 
outbursts, having the same maximum effective temperature at 
the disc edge. Results for model 3 are very similar to those of 
model 2. 

Assuming that L t idai is radiated as a blackbody, the specific 
intensity of the radiation (integrated over frequency) at the disc 
edge is: 



' = vT^Jn 



(3) 



Let a) be the solid angle subtended by the edge of the disc 
at the L\ point; neglecting limb darkening and since the inten- 
sity Eq. [3] is isotropic, the L\ point receives an incoming flux 
F = oT 4 d (w/?r) from the disc edge. The computation of a> 
is implemented in our code and is performed during the time 
evolution of the disc. The ratio F/F+, where F* = oT 4 is the 
intrinsic stellar flux, is given in Table 2. F is equal to a few 
times the intrinsic stellar flux. 



AslSrnak (2002) emphasized, the assumption that tidal dis- 
sipation is radiated entirely at the edge breaks down when r e d ge 
becomes too different from r sur f, the effective temperature at 
the surface of the disc due to viscous heating. In this case, one 
expects an inward radial heat flux, lowering r e d ge and increas- 
ing r sul f . The effective temperature at the surface of the disc (at 
the outer radius) is shown with a dashed line in Fig. 1 . One can 
see that the "viscous temperature" at the outer radius remains 
much lower than the "tidal temperature". Our determination of 
r e d ge is therefore probably overestimated. 

Since the secondary cannot occult the disc completely, an 
enhancement of the disc edge luminosity sh ould be detectable 
durin g eclipses of high inclination systems. IShafter & Misselt 
(2006) analyzed eclipse profiles of the classical nova V Per. In 
this system, the accretion disc is in a hot steady state and ex- 
tends to its tidal radius. Only model profiles with a hot disc rim 
(^edge ~ 6000 - 10000 K) allowed the authors to find the radial 
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Fig. 1. Results for model 1 (left panel) and for model 2 (right panel). The top curve is the absolute V magnitude of the system 
(the origin of our magnitude scale is arbitrary). The bottom curves show r e( j ge (plain line) and the effective temperature at the 
surface of the disc resulting from viscous dissipation (dashed line). 



temperature profile T oc r~ 3 ^ 4 expected in the steady state case. 
This range of temperature is in agreement with the tidal heat- 
ing at the outer edge. In their model, the authors have for the 
first time taken into account the disc thickness which is crucial 
in determining the temperature of the disc edge. In previous 



studie s of eclipse p rofiles of dwarf novae , see e.g. JBruch et al. 



19961 (OY Car) and lHorne & Cookll 19851 (Z Cha), the disc was 
assumed to be "flat" (negligible thickness), which prevented a 
correct determination of the temperature at the disc edge. 

2.2. Heating by scattering of the accretion flux 

2.2.1 . The dwarf nova case: scattering by a wind 

The observation of UV lines with typical P Cyg profiles dur- 
ing dwarf n ovae ou tbursts implies the presence of mass loss in 
a wind (see IWarnerll2003l) . Under the assumption of spherical 
symmetry of the wind, detailed studies of the line profiles have 
shown that the wind acceleration is slow ; the velocity profile 
can be taken as dPrewll987tlMauche & Ravmondll987l) : 



v(r) = vo + (v M -vo)(l-— f 
r 



(4) 



with P — 6, yielding a slow acceleration. As in 
Mauche & Raymond! d 1987b . we take v = 2 x 10 7 cm s" 1 (we 



will show that this parameter is of no importance in our case) 
and Voo = 5 x 10 8 cm s , which is of the order of the escape 
velocity of the white dwarf. 

Photoionisation models of the wind are consistent with a 
mass loss rate in the range 1CT 10 - 1O~ 9 M yr -1 , which trans 



lates into a few percent of the accretion rate dHoare & Drew 
1993b . 

Due to the presence of matter above the disc, a fraction 
of the accretion luminosity is scattered toward h\ and could 
contribute to its heating. We assume in the following that 5% 
of the accretion rate is lost in a spherically symmetric wind. 
The mass conservation law yields: 



Anr p{r)v{r) - aM. 



ace. max 



(5) 




Fig. 2. Side view of the binary system. The innermost part of 
the wind close to the compact object is shielded from L\ by the 
disc and does not contribute to the scattered flux. 



where M accmax is the accretion rate onto the compact object 
during an outburst (see Table 1), a - 0.05 is the 5% mass loss 
and r is the distance to the primary center of mass. From equa- 
tion|5]we deduce the radial density profile of the wind. 

We consider for simplicity a uniform and isotropic opacity 
k - 0.4 (corresponding to Thompson scattering by free elec- 
trons); the optical depth of the wind is then: 



Ek 



(6) 



where S is the total column density of the wind. The value of 
t for models 1 - 4 is shown in Table [3] r is low and the wind 
is optically thin. Consequently only a very small fraction of 
the accretion luminosity is scattered toward L\ . In order to de- 
termine this fraction, we suppose that no "double scattering" 
occurs and we neglect the radial decrease of the accretion lu- 
minosity due to scattering losses. Scattering is handled as an 
emissivity e equal to: 



e(r) = 



I6n 2 r 2 



p{r)K 



(7) 



where r is the distance to the primary center of mass. The flux 
reflected toward L\ is then: 



Iff'-T 



(8) 
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Table 3. Results for the heating of the secondary by scattering 
of the accretion luminosity, r is the optical depth of the outflow 
and F is the heating flux impinging on the L\ point. 



Model 



F/(L m 



0.01 

0.1 
0.3 



KAnD 1 )) 
3 



F/F+ 



6x 10 
6 x 10- 5 
8 x 10~ 5 



2.1 x 10 - 
2.1 x 10~ 2 
2.6 x 10~ 2 



where d is the distance between L\ and the running point r. We 
show in the appendix that: 



I 



e(r)g(r)dr 



(9) 



where g(r) is a geometrical weight function for the spherical 
shell at radius r. Note that the radial integration starts at r c (see 
Fig. [2]), and not at r\. This approximately takes into account 
the screening of the inner region of the wind by the edge of 
the accretion disc. r c is computed by considering that the disc 
has its maximal radius (see previous section) and typical h/r 
ratio at the edge, see Table [TJ Note that due to this screening, 
the value vo in Eq. [4] which controls the mass density of the 
wind near the white dwarf, does not influence the fraction that 
is back-scattered. However, it does influence the value of E and 
hence r, so the values of r given in Table [3] should be taken 
with care. 

Equation [9] is numerically computed with the appropriate 
e(r) profile. In the case of a direct illumination, the L\ point 
would receive a flux equal to L acc /4^D 2 . As shown in Table 3, 
the fraction of this flux scattered toward L\ is very small, of the 
order of 1CT 4 . This yields a very low F/F* ratio. DNe have a 
insufficient luminosity for this effect to be significant. 

2.2.2. The X-ray transient case: scattering by a corona 

Observational evidence for an accretion disc corona in 
low mass X- r ay bin arie s (LMXBs) was f i rst di scussed by 
White & Holt] ( 1 1982b and iMcClintock et alj dl982b . The prin- 
cipal indicator of an extended region of hot, ionized gas above 
the disc is the evidence of partial eclipses of X-ray emission in 
high inclination systems. While the compact object is screened 
from the observer, a fraction as high as 50% of the X-ray flux 
is left during eclipses. These sources are the so-called "coronal 
sources". Such a corona is believed to be present in many, if 
not all, LMXBs. X-ray transients in outbursts should not be an 
exception. 

We discuss here the possibility that, during an outburst, the 
L\ point receives some fraction of the accretion flux due to 
back-scattering by the corona. 

The geometry of the accretion disc corona is still poorly 
understood. One can grossly distinguish between two types of 
geometry: a slab corona which consists of two slabs of opti- 
cally thin plasma sandwiching the accretion disc or a spherical 
corona located around the compact object. 

Let us consider first a spherical corona. With typical X-ray 
transient parameters (see TableQ}, the value of r c (see figure[2| 
is r c ~ 10 11 cm. According to studies of eclipse light-curves 



of coronal sources, the most extended coronae are found in 
the brightest sourc es (Ly ~ 10 39 erg s' 1 ) and have a radius 



of 7 x 10 1() cm (see lChurch & Balucihska-Churchl2004l) . Such 
a corona can be screened by the accretion disc and would not 
scatter radiation toward L\ . 

A slab corona results from the X-ray heating of the accre- 
tion disc surface by the central source. The heated gas forms 
a tenuous layer with a thickness exceeding that of the disc. 



Begelman et al.1 (11983b analyzed the structure and dynamics of 



the X-ray heated disc. They found that, in addition to the for- 
mation of the corona, X-ray heating could also drive a thermal 
wind in the outer region of the disc. A slab corona could back- 
scatter radiation coming from the central part of the disc over 
the disc rim, toward the L\ point. A wind emitted from the outer 
region of the disc would also contribute to this effect. It is how- 
ever hard to obtain a quantitative measure of this effect. 

Results obtained in the study of irradiation of the accretion 
disc in X-ray transients can give us a clue about the efficiency. 
iDubus et al.l 0999) show that self-consistant computation of 
the radial structure of irradiated discs in SXTs leads to a prob- 
lem: the disc has a convex shape and the outer, unstable, regions 
are screened from irradiation. This is contrary to observation 
and the authors replace the usu al formula for the irradia tion 
temperature in the disc (see e.g. lShakura & Svunvaevl 1973b by 
the following prescription: 



art 



c 



Anr 2 



(10) 



where C is taken to be constant, which is thought to reflect 
the fact that the disc receives some X-ray flux at each radius. 
This could be due to warping of the disc, or to the presence 
of a corona. In SXT, C — 5 x 10" 3 is an adequate value (see 
Dubus et al.lll999l EOOlb . Note that a warping of the disc in 
its outer regions could allow a direct irradiation of L\ . If the 
corona is responsible for the overall heating of the disc, C could 
serve as a first estimate for the efficiency of the heating of the 
L\ point by the corona. The flux impinging the L\ region would 
therefore be: 



F = C 



4nD 2 



(11) 



where D is the distance between Ly and the compact object. 
SXT outbursts have typical luminosities that translate into an 
irradiation flux at the secondary of the order of 10 5 F*. With 
C — 5 x 10~ 3 , Eq. QTJyields F - 500.F* which is very signif- 
icant. This value is most probably overestimated ; a better de- 
termination, relying on a model for the corona and/or the wind, 
should be sought. 



3. Variation of M tr 

3.1. Vertical structure in quiescence 

The mass transfer rate is given by (see lLubow & Shulll975h : 



Mr = Qp(L x )c s 



(12) 
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where c s is the isothermal speed of sound which depends on 
the temperature at L\, T(Li): 



c s = ^8,r(Li) (13) 

and Q is the cross section of the stream (ILubow & Shull 1975b : 



^ 



e- 2 iW 



(14)1 



where k depends only slightly on the mass ratio q. Here we take 
k-1, valid for q in the range 0.1 - 0.6 within 10 %. 

In order to determine p(L\) when an incident flux heats the 
upper layer of the atmosphere, we first compute the structure 
of the envelope of the secondary during quiescence for models 
1 - 4 by solving the standard equations: 



10000 r 



000 




10 z 10° 10 4 10 5 



1 D 1 7 1 0° 1 s 
P (dynes. cm -2 ) 



dp ^ 

d\nT 

= V 

dlnP 

P = P(P,T) 



(15) 

(16) 
(17) 



where V is either the radiative temperature gradient when V ra< j 
is less than the adiabatic value V a( j with 



V rad = 



kPF 

4PradCg 



where P la( j is the radiative pressure, or the convective value 
when the radiative gradient is superadiabatic (V ra( j > V a( j), 
calculated in the mixing-length prescription as described in 
Paczvriskil dl969l) . The equat ion of state and opacities are com- 
puted in the same way as in lHameurvl d 1991D . We use the diffu- 
sion approximation for the radiative flux, which is not accurate 
in optically thin regions. The true radiative flux should be com- 
puted by solving the radiative transfer, a task which is beyond 
the scope of the present exploratory work. The photosphere of 
the atmosphere is defined to be the location where r = 2/3, 
where r is the Rossland optical depth. 

We solve numerically Eqs. ( [T5lfT7] i assuming a given T(L{) 
and p(L\) at the L\ point. T(Li) is chosen so that the tempera- 
ture of the model at the photosphere is equal to the value of T* 
given in Table 1. The value of p{L\) is then determined from 
Eq. (fT2b with the appropriate value of T(L\) and with the secu- 
lar mean value of M tr given in Table 1 . The computation starts 
at the L\ point and goes inward until a temperature > 10 4 K 
is reached. The Roche geometry is taken into account by using 
the Roche potential in the hydrostatic equilibrium equation. 

Figure 3 shows the vertical structure of the secondary en- 
velope for DNe below the period gap (model 1), DNe above 
the period gap (model 2 & 3) and SXTs (model 4). As seen in 
Fig. |3] the L\ point lies in the radiative part of the atmosphere 
in the SXT case and in the convective part of the atmosphere 
in DNe secondaries. In each case, the L\ point is located above 
the photosphere, with zl, /H ~ 2 in model 4, Zl, /H < 1 in 
model 1&2 and zl x /H ~ 0.3 in model 3 (H is the vertical pres- 
sure scale height of the envelope). This difference is due to the 
fact that SXTs have lower mass transfer rates than DNe, which 
is possible since the critical mass transfer rate above which the 



Fig. 3. Vertical structure of the envelope for a secondary in DNe 
below the period gap (green curve) and above the period gap 
(blue curves) and for a SXT (black curve). All curves start at 
the L\ point and end at an arbitrary depth. The line is dashed 
when energy transport is radiative. The dots mark the position 
of the photosphere, defined to be the location where r = 2/3. 
The models are constructed in order to have T(j - 2/3) = T*. 



(18) instability disappears i s lower in an irra diated disc than in the 



unirradiated case (see lDubus et al.lll999l) . 



Figure|3]also shows that for DNe below the period gap, the 
temperature increases very slowly with depth in the secondary 
envelope. This is due to the low value of the adiabatic gradient 
(V a( j ~ 0.1) in the regions of molecular hydrogen dissociation 
present in the envelope of these low temperature stars. For DNe 
above the period gap, the secondary is hotter and the L\ point 
lies in a region where hydrogen is neutral. As a consequence, 
the adiabatic gradient is larger (V a( j = 0.4) and the temperature 
increases rapidly with depth. 

3.2. Vertical structure during an outburst 

We now assume that an heating flux Fh ea t strikes the L\ region 
and is absorbed below the photosphere. i*h ea t heats up the sec- 
ondary on the thermal time scale of the atmosphere, short com- 
pared to the duration of outbursts, until the envelope reaches an 
outward flux equal to: 



o-Tl + /3F } 



heat 



(19) 



Here T* is the effective temperature of the secondary in quies- 
cence and (3 is the bolometric albedo that accounts for a hor- 
izontal transport of the deposited energy, which could be im- 
portant in the case of deep convective stars (for radiative stars 



f5 — 1, see e.g. ISmithll 19951) . Since its value is not known p re- 



cisely, we adopt/? = 1 in all cases (see lBrett & Smith 1993). It 
is not clear whether on the time scale of an outburst an efficient 
horizontal heat transport can set in. 

The evol ution of the sub photospheric layers has been de- 
termined by lHameurv et al.l (1988). As energy is deposited in 
sub-photospheric regions, heat diffuses inward and hinders the 
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Fig. 4. Temperature profiles in models of atmospheric heating by 100.F+. Model 3 yields the same results than model 2. 



escape of the intrinsic stellar flux. An isothermal layer of tem- 
perature ri ayel = (F out /cr) 1 ^ 4 form, on top of matter that has 
not yet been heated. This heated layer ultimately extends to the 
depth where the quiescent temperature profile is equal to ri ayer . 
However, assuming that the layer has been heated entirely by 
the stellar flux, its depth is restricted by the requirement that: 



Eth < Afou^m-stC-^* 



(20) 



where At is the duration of the outburst and E± is the thermal 
content of the layer; we take here At = 30 P OI b for models 1-3 
and At = 100 P olb for model 4. We then solve Eqs. dl54417T > 
to determine the density profile of the isothermal layer, and the 
new mass transfer rate is then found from (Eq.fT2l with the new 
values of T(L\) andp(Li). 
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Fig. 5. Mass transfer enhancement in SXT (model 4, black 
line), DN above the period gap (model 2 & 3, blue line) and 
DN below the period gap (model 1, green line). 

Examples of models where the heating flux is lOOF* are 
shown in figure |4] The figure shows that in this case, the 
isothermal layer extends down to a column density of the or- 
der of 10 2 g cirT 2 in dwarf novae and 10 3 g cirT 2 in soft X-ray 
transients. The depth of the layer is lower in dwarf novae be- 
cause of the restriction (|20T > in low period systems (SU UMa) 



and because of the location where T - T\ aya is quickly reached 
in longer period systems (U Gem and Z Cam). 

The mass transfer enhancement is computed for a heating 
flux Fheat ranging from 0.1F* to 100F+. The results are shown 
in Fig. [5] As can be seen, the mass transfer rate enhancement is 
quite significant in SXTs (factor of 2 for Fheat = 0.1F*) and is 
potentially very significant (factor of 100 for Fh ea t = 10F*). 
In comparison, the mass transfer enhancement is moderate 
in DNe. This is a direct consequence of the different depths 
reached by the isothermal layer. 



3.3. Absorption of the incident radiation by the 
secondary atmosphere 

The effect of an incident flux on the secondary atmosphere de- 
pends on the spectrum of the incoming radiation, and hence 
on its origin. The flux absorbed below the photosphere, de- 
noted Fheat in the last subsection, can be quite different from 
the incident flux F. In the case of heating of L\ by the hot 
edge of the disc, the spectrum of the impinging radiation can be 
aproximated by a blackbody spectrum with effective tempera- 
ture =* 10 4 K. In the case of heating by scattering of the accre- 
tion luminosity, we assume for simplicity that no reprocessing 
occurs. As one half of the accretion luminosity is emitted in 
the accretion disc and the other half is emitted in the boundary 
layer, we use two separate contributions to the energy spec- 
trum of the incident radiation. The boundary layer emission is 
assumed to be a black-body ; Table [4] shows the wavelength of 
the emission peak for each model. For the disc emission, we 
use a standard multi black-body spectrum of a steady disc, in 
which we have taken into account the irradiation by the central 
source in the case of SXTs. Steady state is a very reasonable 
assumption for discs during outburst. 

Radiation absorbed in optically thin layers contributes to 
the formation of a hot corona above the photosphere ; only the 
incoming radiation that penetrates below the photosphere con- 
tributes to changing the vertical structure and hence increasing 
the mass transfer rate. We use the ATLAS 12 code to com- 
pute the monochromatic opacities in the envelope and deter- 
mine which fraction reaches the photosphere. The ATLAS 12 
opacities extend from 10 /mi down to 10 nm, which is a suitable 
range for our purpose except for the SXT case where the bulk 
of the accretion flux is emitted in the soft X-ray band (emission 
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Table 4. Fraction of the incident radiation emitted by the disc edge, by the accretion disc and by the boundary layer reaching the 
photospheric level (t = 2/3). See text for details. 



Model 


Disc edge at T eff = 10 4 K 


Disc 




Boundary layer 


1 

2,3 
4 


0.25 
0.18 
0.22 


0.01 
0.05 
0.1 


A pic (nm) 

50 

30 

1 


quiescence /Vat/^* = 100 
10 2 0.1 
< 10~ 3 2 10~ 2 
~0 0.5 



peak at ~ 1 nm). In this last case, the opacity is extrapolated 
down to very short wavelengths and includes Thomson scat- 
tering. Both continuous and line opacities of ATLAS 12 are 
described in Castelli (2005). Since we consider here stars with 
effective temperatures below 5000 K, molecular opacities are 
also included. 

Figure [6] shows the monochromatic opacity at the photo- 
sphere of model 2. The opacity profiles of the other models 
show qualitatively the same behavior, the most important com- 
mon feature being the very high opacity below -100 nm due 
to neutral hydrogen. 

We define the monochromatic optical depth t v : 



Tvfe) 



I 



iLO 



pK v dz 



(21) 



We start at L\ with an incident specific intensity B Y (T) (black- 
body of temperature T) and we compute the specific intensity 
at each depth z that results from absorption: 



I v (z) = e- TAz) B v (T) 



(22) 



Finally at the photosphere (depth z p ) we compute the frac- 
tion of energy left: 



a( TR = 2/3) = 



J Iy(z p )dv J Iy(z p )dv 
JB v (T)dv ~ o-T 4 / n 



(23) 



The results are shown in Table|4] For the radiation coming from 
the disc edge, approximately 20 - 25 % of the flux is able to 
penetrate below the photosphere. 

The radiation emitted by the boundary layer is almost 
completely absorbed above the photosphere of the secondary. 
This is because the black-body maximum is located below the 
Lyman limit (A ~ 90 nm) and absorption by neutral hydr ogen 
is ver y important. This point was already discussed by King 
dl989 ) in the context of cataclysmic variables. Between 1 % and 
10% of the radiation emitted by the accretion disc reaches the 
photosphere. This is because the outer region of the accretion 
disc emits at lower energy than the boundary layer, i.e. above 
90 nm. One can check that the flux left at the photosphere is 
mainly in the optical and near infra-red bands. 

Note however that our computation here is oversimplified 
since we do not consider at all the alteration of the atmosphere 
by the incoming radiation. 

We have also computed a using heated profiles we con- 
structed in Sect 3.2. Our results show that the column den- 
sity of the photosphere decreases with an increasing heating. 




100 



1000 



10000 



A (nm) 



Fig. 6. ATLAS 12 opacities at the photosphere of model 2, 
where p ~ 10~ 7 g cm~ 2 and T - 3500 K. The line profiles 
have been smoothed for readability. The dashed line is the con- 
tinuum opacity. 



As a consequence, an increasing fraction can be absorbed be- 
low the photosphere. Table |4] shows the fraction of the inci- 
dent flux absorbed in the unperturbed atmosphere of the sec- 
ondary, and that absorbed in an atmosphere efficiently heated 
with Fheat/^* = 100. We do not discuss changes for dwarf no- 
vae, since the scattered flux is far too small to have any effect. 
However, model 4 shows a huge difference: in quiescence noth- 
ing reaches the photospheric level whereas in the heated state, 
the photospheric level is so close to L\ that most of the incom- 
ing flux reaches the photosphere. We therefore cannot deduce 
from such simple arguments whether the heating could be effi- 
cient or not. To obtain a quantitative result, one needs to solve 
the full radiative transfer problem. This is left for a future in- 
vestigation. 

4. Conclusion 

We have considered two mechanisms that could be responsi- 
ble for heating of the L\ point during outbursts of dwarf no- 
vae and soft X-ray transients. First, the outer edge of the disc 
is heated by tidal dissipation and reaches an effective temper- 
ature of the order of 10 4 K (see Sect. 1). We have shown in 
Sect. 3 that a significant part of the resulting thermal radia- 
tion can penetrate below the photosphere of the secondary. The 
disc edge is able to heat the upper layer of the secondary with 
an incoming flux comparable to the intrinsic stellar flux. This 
has probably no significant effect in dwarf novae. In soft X-ray 
transients, the dominant effect could be a heating by scattered 
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radiation. A disc corona, or a wind driven by the X-ray heating 
of the disc, could act as a scattering medium, allowing the radi- 
ation to overcome the screening by the accretion disc. It is hard 
however to infer a quantitative estimate of this effect. However, 
the luminosity of SXTs are so high that even a low efficiency 
(1CT 3 - 1CT 4 ) could lead to a significant heating. However, the 
bulk of this radiation is absorbed by neutral hydrogen, and ini- 
tially nothing reaches the photospheric level of the secondary 
star. Nevertheless, the flux is so strong that the incident radi- 
ation probably strongly modifies the state of the atmosphere. 
Yet it is hard to assess whether in this case an efficient heating 
takes place or not. The resolution of the corresponding radia- 
tive transfer should be undertaken. 

In Sect. 3 we determined the relation between the mass 
transfer enhancement and a given heating flux impinging the 
L\ point (independently of the origin of this heating). We com- 
puted the vertical structure of the secondary with the assump- 
tion of hydrostatic equilibrium. The atmosphere is not hydro- 
static near the L\ point due to the mass outflow in this region. 
However, departure from hydrostatic equilibrium is likely to 
be important only in the near vicinity of the L\ point. We found 
that the L\ point lies higher in the atmosphere in soft X-ray 
transients than in dwarf novae. This is a consequence of a 
lower value of the mass transfer rate in soft X-ray transients, 
where the L\ point lies in a nearly isothermal, radiative region 
of the atmosphere, whereas in dwarf novae the L\ point lies 
in the convective part of the atmosphere, where temperature in- 
creases more rapidly with depth. As a consequence, an incident 
flux can heat a significant layer in soft X-ray transients, which 
yields a large enhancement of the mass transfer rate. On the 
other hand, in dwarf novae only a narrower region of the atmo- 
sphere is affected. For a dwarf nova below the period gap, the 
envelope is too massive to be fully affected on the time scale 
of an outburst. Our results suggest that the mass transfer en- 
hancement could be potentially very important for soft X-ray 
transients, up to a factor of ~ 100 for a heating flux equal to 
10F*. For dwarf novae, the mass transfer rate enhancement is 
more moderate, of a factor of ~ 10 for a heating flux equal 
to lOF*. In dwarf novae, none of the effects investigated here 
is able to produce such fluxes. The soft X-ray transient case 
needs further investigation, on one hand to determine the frac- 
tion of radiation that is back-scattered toward L\, and on the 
other hand the efficiency of the heating of the secondary has to 
be determined properly. 
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Appendix A: Derivation of Eq.|9j 

We give here the detailed computation of the integral: 
e(r)dV 



fir 



d 2 

with 



(A.l) 



e(r) = J6& p(r) * and P^^P^f' 2 ^ A - 2 > 

We adopt a spherical system of coordinates centered on the 
primary center of mass, with the colatitude coordinate 6 de- 
fined to be the angle between the running point r and the line 
joining the primary center of mass to L\ and with longitude co- 
ordinate (p. Due to the symmetry of the problem, the integral 
lA.llbecomes: 



£ dMr) ff 



dS 
If 



dre(r)g(r) (A.3) 



with dS = r 2 sin 9d6d<p and d 2 - r 2 + D 2 - 2rD cos 0. Note 
that the radial integration starts at r - r c (see Fig. [2) and not 
r — r\ this approximately accounts for the shielding of the scat- 
tering region by the disc rim. The double integral involving the 
angular variables gives: 



rCdS r 1 + 7 

g (r) = — = 2tt- In r - (A.4) 

* JJ d 2 D |l - 2| 

It is useful to look at the asymptotic behavior of this func- 
tion: one can check that g(r) — » 47r as r — > oo and that g(r) — > 
as r — » 0. Note that this quantity diverges at r — D, however 
the integral Eq. |A.3| is well defined because the divergence of 
the logarithm function at the origin is integrable. 



